Interaction-driven transition between topological states in a Kondo insulator 
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Heavy fermion materials naturally combine strong spin-orbit interactions and electronic correla- 
tions. When there is precisely one conduction electron per impurity spin, the coherent heavy fermion 
state is insulating. This Kondo insulating state has recently been argued to belong to the class of 
quantum spin Hall states. Motivated by this conjecture and a very recent experimental realization 
of this state, we investigate a model for Kondo insulators with spin-orbit coupling. Using DMFT, 
we observe an interaction-driven transition between two distinct topological states, indicated by a 
closing of the bulk gap and a simultaneous change of the Z2 topological invariant. At large inter- 
action strength we find a topological heavy fermion state, characterized by strongly renormalized 
heavy bulk bands, hosting a pair of zero-energy edge modes. The model allows a detailed under- 
standing of the temperature dependence of the single particle spectral function and in particular 
the energy scales at which one observes the appearance of edge states within the bulk gap. 

PACS numbers: 71.27.+a, 71.10.Fd, 03.65.Vf 
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I. INTRODUCTION 

Following the theoretical discovery of the Quantum 
Spin Hall (QSH) effect^"— and its experimental realization 
in HgTc quantum well systems^ the field of topological 
insulators has attracted a lot of research interest. Most 
intriguing is the fact that this state cannot be adiabati- 
cally connected to a conventional insulating state with- 
out breaking time reversal symmetry (TRS). As a con- 
sequence, on the boundary between a topological insula- 
tor and a trivial insulator a pair of helical edge modes, 
which are robust against disorder and interaction effects, 
emerges^. The interplay of correlation effects and topo- 
logy has become a very active field of research^. Recently 
a study of the so-called Kane- Mele- Hubbard models by 
means of quantum Monte Carlo revealed, that the topo- 
logical insulator phase exists up to rather large values 
of U before undergoing a continuous phase transition to 
an antifcrromagnctically ordered state, thereby breaking 
TRS. In the disordered phase the interplay of the collec- 
tive spin mode and topology of the band structure leads 
to novel effects and model system s 10 ' 11 . 

Hubbard-like correlations have equally been included 
within dynamical mean-field theory (DMFT) for a va- 
riety of models including the BHZ model^ of the QSH 
effect in HgTe quantum wellsH^. At the DMFT level, 
collective modes are absent and correlation-driven tran- 
sitions can be understood in terms of a renormalization 
of the band structure due to a constant real part of the 
self energy or due to dynamical effects, which can lead to 
the divergence of the effective mass. 

DMFT is the method of choice to understand the 
salient temperature dependent features of the paramag- 
netic heavy fermion stated. Below the coherence tem- 
perature, the individual Kondo screening clouds of the 
magnetic impurities overlap coherently to form the heavy 
fermion liquid. This state of matter is adiabatically con- 
nected to the non-interacting system and is hence a Fermi 



liquid. The ultimate signature of coherence in models of 
heavy fermions is the Kondo insulating state, where there 
is precisely one conduction electron paired with an im- 
purity spinii. Prototypical Kondo insulators are YbB 1 2 , 
CeNiSn, Ce 3 Bi 4 Pt 3 and SmBe^. As proposed in Ref. |l| 
and a state with a non-trivial band topology may be 
realized in Kondo insulating systems due to the strong 
spin-orbit coupling present in these materials, even in the 
presence of strong electron-electron correlations. In other 
words, under certain conditions the Kondo insulating 
state is adiabatically connected to a non-interacting insu- 
lator with a non-trivial band topology. This results from 
the odd-parity wave-function of the f-electrons^ -, which 
leads to a non-trivial, momentum-dependent hybridiza- 
tion between f-electrons and conduction electrons. Very 
recently, the Kondo insulator SmBg was investigated us- 
ing an experimental setup suited to identify non-local 
transport effects^. A topological Kondo insulating state 
was indeed found to exis1i2i~— , explaining very naturally 
both the residual conductivity as T — > and the exis- 
tence of in-gap states found in ARPES measurements^. 
Follow-up investigations using a similar method may very 
well find topologically non-trivial states in other Kondo 
insulating materials, which are up to now not very well 
understood^!. 

The aim of this paper is to investigate models of heavy 
fermions within the DMFT framework which have a topo- 
logical Kondo insulating ground state. One of our central 
interests is to assess if this state of matter can occur in 
the strong coupling local moment regime. In section |TT] 
we describe our model for topological Kondo insulators. 
Next we briefly introduce the concept of the topologi- 
cal invariant and explain its evaluation for an interacting 
system, using the so-called topological Hamiltonian. Af- 
ter describing our method in section IIV1 we determine 
the topological properties of the model. In sections IVl 
and I VII we present our main findings, before concluding 
in section IVfTl 
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II. TOPOLOGICAL KONDO INSULATOR 



the QSH effect. It is given by 



We will devise a model appropriate for Cc-based com- 
pounds in 2D2£. Our starting point is the Ce ion's J=5/2 
multiplet, which originates from the 4f-orbitals in the 
presence of spin-orbit coupling. In the solid state, the 
multiplet is split further by the crystal electric field, with 
the splitting depending on the lattice symmctr y22. We 
assume that the ground state of the magnetic ion is a 
Kramer's doublet, which we call f-states. In the sim- 
plest case, they hybridize with a single conduction band 
(d-states). With these ingredients, we build a model 
quite similar to the Periodic Anderson Model (PAM) on 
a 2D square lattice. The resulting Hamiltonian is 



H{k) 
T-(-BHz(k) 



(U B Hz{h) 

o n* BHZ 



■(-*>)) 

— cos(fc x ) — cos(fcj,)^e 
- \(sin(k x )cr x + sm(k y )(Ty 



(2) 



which can be obtained from ([1} in the mixed- valence limit 
tf = —t. The parameters are related by A = V and 
m = —Ef/A. Thus we expect some similarities especially 
in the weak coupling limit. Our primary interest however 
is the correlation-dominated Kondo insulating state. 



(1) 
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where <&(k), Ed(k) and Ef(k) are 2x2 sub-matrices, and 
encode all the information on the geometry and the ef- 
fect of spin-orbit coupling. We use (d\ fl) as a short 

hand notation for (dj^. /jtj_), creating an electron 

in the conduction band and the almost localized band, 
respectively. Here t" and 4- denote the pseudo-spin quan- 
tum number. 

Due to spin-orbit coupling the f-states are eigenstates 
of the total angular momentum J, and hence hybridize 
with conduction band states with the same symmetry. 
This gives rise to the momentum-dependence and non- 
trivial orbital structure of the form-factor <&(fc)2£. Fol- 
lowing the derivation in Ref. [26l . we obtain 

Ed(k) — — 2i(cos(fc x ) + cos(fc y ))l 

E f (k) = (sf -2tf(cos(k x ) + cos(k y )))l 

$(&) = d(k)oa 
where d(k) = (2 sin(k x ), 2 sin(A^), 0) 

formulated in the basis of the Pauli matrices <x and the 
2x2 unit matrix 1. The fact that Edj(k) (d(k)) is even 
(odd) under the transformation k — > —k guarantees time 
reversal symmetry. The interaction term Hu acts only on 
the almost localized f-states. The delocalized conduction 
electrons are only weakly affected by correlation effects, 
therefore we neglect interactions on the d-orbitals. We 
let the hybridization amplitude take a moderate value 
of V = 0.4 1, and allow for a small, hole-like hopping 
amplitude between f-orbitals, i.e. tf = —0.2t. Both will 
be renormalized to a smaller value in the presence of the 
interaction. 

The non-interacting part of our model can be related to 
the BHZ mode l 12 ' 14 for HgTc quantum wells exhibiting 



III. TOPOLOGICAL INVARIANT 

The topological insulator docs not connect adiabati- 
cally to the trivial band insulator—. Both states cannot 
be characterized by a local order parameter and one has 
to resort to the concept of a global Z2 topological in- 
variant v = 0, 1 to distinguish theroi. Provided that the 
protecting time reversal symmetry is not broken, states 
characterized by different values of v can't be adiabati- 
cally connected without closing the single particle gap. 
A formulation of the topological invariant in terms of 
the single-particle Green functional quite naturally ex- 
tends the concept to interacting systems adiabatically 
connected to the non-interacting case 2 ^. It turns out 
that the zero-frequency value of the Green function suf- 
fices to characterize the state in terms of the topological 
invariant 2 ^, which leads to the idea of defining a so-called 
topological Hamiltonian^ 

hto P o(k) = -G -1 (Ma; = 0) = h (k) +E(/c, iw = 0). (3) 

The topological invariant of the interacting system is 
given by the topological invariant corresponding to 
the non-interacting system governed by the topological 
Hamiltonian h topo . It characterizes correctly the system 
in terms of the topological invariant, but it doesn't pro- 
vide realistic quasi-particle spectral. In the presence of 
inversion symmetry, the calculation of the topological in- 
variant of §3§ is even further simplified^. One simply 
has to calculate the quantity 



(4) 



where m is the band index, Ti are the four time reversal 
invariant momenta in 2D, and (2™ is the corresponding 
parity eigenvalue at 1^ of the 2m-th band. The Zi topo- 
logical invariant is then obtained via 



(-ir=n* 



(5) 
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IV. METHOD 



formulated in terms of an effective action 



We employ the DMFT to make our model system ([T]) 
accessible to numerical methods 3 ^. The central approxi- 
mation of the DMFT is the assumption of a momentum- 
independent self energy 



E(fc, iuj) ~ E(iu>). 



(6) 



This approach neglects non-local correlations. In addi- 
tion, as we are interested in the properties of the para- 
magnetic phase, we do not allow for magnetically ordered 
states, which would break TRS. This translates to a sim- 
plification of the topological Hamiltonian ([3]), which is 
now the original Hamiltonian renormalized by a constant 

hto P o(k) = ho(k) + S(»w = 0) = h (k) + E . (7) 

We can now find an explicit formula for the Zi topological 
invariant ((5|) in terms of our model's bare parameters and 
the renormalization constant Eo: 

(— 1)" =sign(e / + E ) 2 x sign(4t - % - E -£/)x 
ign(4i/ -4i-E -£/) 

(8) 



x Slf 



The resulting phase diagram in terms of £/ and Eo is 
shown in Fig. [TJ It represents a roadmap, which sim- 
plifies the characterization of the system's state once 
Eo = Eo(£/) is known. The simplification allows for a 
mapping of the lattice problem (QJ to an auxiliary single- 
impurity problem, which describes a fully correlated im- 
purity Wimp coupled to a non-interacting bath. It can be 
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FIG. 1: Phase diagram of the topological Hamiltonian ((T)) 
for tf = — 0.2 1. On the lines where the right-hand side of 
equation (|8]) evaluates to zero, the topological invariant is 
not defined, and the system is in a gap-less, semi-metallic 
state. The colored, thick lines represent the simulation runs 
for different Ef and U and the analytical results at U = 0, 
were Eo = 0. The line color corresponds to the band gap 
size A 9 . The actual gap closings coincide very well with the 
predicted transition lines. 



S eff 



J dTdT'p(r)A(T-T / )f(r') + J drU lmv {r). 



(9) 

The bath parameters A(r) have to be obtained in a self- 
consistent manner. Once convergence is reached, the im- 
purity self-energy is equivalent to the local self-energy 
E(iw) of the original system. 

We use the numerically exact CT-HYB quantum Monte 
Carlo algorith m 33 ' 34 as the impurity solver at low 
but finite temperature, which yields the self-energy in 
imaginary frequency. The inverse temperature was cho- 
sen as fit = 100 in most cases, while close to the transi- 
tions it was increased up to fit = 300. 
In the Fermi liquid regime, Im[E(iw)] ~ iuj, while the 
real part can be extrapolated reliably to a finite value, 
i.e. Re[E(iw — > 0)] = Eo- From the local one-particle 
Green's function in imaginary time G/ oc (t) ~ exp(re±) 
we infer the position of the upper and lower band edges 
e±. 



V. TRANSITION BETWEEN DISTINCT 
TOPOLOGICAL STATES 

In Fig. [1] the phase diagram of the topological 
Hamiltonian is shown, superimposed with numerical 
results for the evolution of the band gap size A g as a 
function of Eo(J7), for different values of Sf. Positioning 
the bare f-level at Ef = — 2.0 i puts the system in a 
topological insulating state at U = 0, with a bulk band 
gap A g sa 0.88 1. On a geometry with open boundaries or 
with an interface to a conventional insulator, this leads 
to the formation of a pair of zero-energy edge modes 
inside the bulk gap. They cross at the T-point k x = 0, 
as can be seen in Fig. [2^. We call this state the (0)-TI 
state. As we slowly increase U, the system stays in this 
phase as long as the bulk gap remains open and TRS 
is not broken spontaneously. The results for the bulk 
band gap size A g arc shown in Fig. [3] While increasing 
U, we adapt the position of the chemical potential fi to 
always stay approximately in the center of the band gap. 
This way we also ensure, that the system is precisely 
half-filled. Approaching U = 4.0 1, the size of the band 
gap smoothly decreases, and the gap eventually closes 
at U = 4.0 1. Right at this point the system is in a 
semi-metallic state where the topological invariant is not 
defined, and can therefore change discontinuously while 
crossing this point in phase space. In the phase diagram 
of Fig. Q]this point is exactly on the black, singular line 
given by Eo = — £/. Increasing U even further leads to 
a reopening of the band gap. The calculation of the 
topological invariant using the topological Hamiltonian 
([7]) reveals, that for U > 4.0 1 the system is still in 
a topological insulating state. Not surprisingly, it is 
distinct from the state at [7 = 0. Due to the band 
gap closing at U = 4.0 1, an adiabatic connection to 
the (0)-TI state is not possible. Instead this state is 
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FIG. 2: Band dispersion for the non-interacting case with 
open boundaries in y-direction, N y — 24, a) crossing at F, b) 
crossing at X. 




FIG. 3: Evolution of the bulk band gap with increasing U for 
different values of £/. The solid curves are guides to the eye. 



connected a different topological state, namely the one 
present for £/ > 0. We call this state the (7r)-TI state, 
as in this case the edge states cross at k x = ±n, as 
shown in Fig. ^jp. The adiabatic connection can indeed 
be established by smoothly changing e / and U in a way, 
such that the system stays in this phase. This is done 
in a simulation run starting at Sf — — 2.0 i, U = 6.0 1, 
which corresponds to Eo ~ 2.3 1, and evolving the 
parameters to £/ = +1.0 1 and U = 0, while staying 
well away from the singular line Eo = — £/. Along 
this path we find no gap closing, as the upper-most 
colored, somewhat wiggly line in Fig. [T] clearly shows. 
This signifies that at U = 4.0 1 we observe a transition 
between two distinct topological states, which is driven 
by the interaction. 



VI. HEAVY BAND TOPOLOGICAL 
INSULATOR 

One can think of other paths in the phase diagram of 
the topological Hamiltonian. For the choice £/ = — 6.0 1 
the non-interacting system is an insulator with a trivial 
band topology. The weakly dispersing f-band is separa- 
ted from the d-band by a band gap A g ps 1.2 1. Smoothly 
increasing U leads to a phase transition to the (0)-TI 
state once the transition line in Fig. [1] defined by 
Eo = 4tf — At — £f is reached, which is the case at 
U ps 1.25 i. Indeed, this transition is the one observed 
in Ref. and [H in the BHZ-Hubbard model for m > 2 
at weak coupling strengths. Beyond the semi-metallic 
transition point, as the system is driven deeper into the 
topological insulator phase, the band gap size increases. 
It starts to decrease again beyond U = 3.0 1. For values 
of U > 8.0 1 the bulk band gap size becomes very small, 
signaling the proximity to the next transition line, which 
separates the two distinct topological insulator states of 
(0)-TI and (7r)-TI. Fig. 2] reveals that as we increase U 
beyond the value of U = 8.0 1, the evolution of the renor- 
malization constant Eq(J7) becomes very flat. While it 
approaches the value of Eo = 6.0 1, where the transition 
should happen, the effective mass, as defined by 



= 1 



<iImE(u;) 



dui 



wi _LnE^n) 



increases rapidly, and apparently diverges at U ps 8.9 1. 
As a consequence, Im[E(zo;)] remains finite in the limit 
uj — > 0, indicating a possible orbital-selective Mott 
transition of the heavy band due to local dynamical 
fluctuations^, or a transition to an antiferromagnetic 
state, which can't be captured in our approach. The 
possibility of a TRS breaking transition is supported by 
the fact that the f-band is exactly half- filled at this point, 
while the double occupancy is strongly suppressed. The 
nature of the transition and the phase lying beyond the 
transition can't be investigate with our current method. 
Most probably, due to the approximation of DMFT, in 
our simulations the size of the band gap becomes so small 
that it can't be resolved, but it never close o 14 ' 37 . 

To quantify the importance of correlations, and espe- 
cially to pin down the local moment regime, it is inter- 
esting to investigate the quantity 
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(n 



(/)„(/) 
t n l 



(n\ f) )(n[ f) ) 



111 



where = without correlations. For a local moment, 
where the double occupancy is completely suppressed, 
= 1. Thus we can use this quantity as a measure 
for the impact of correlations on the system. The evo- 
lution of function of U is plotted in Fig. [5] 
With increasing U, the occupancy of the f-states de- 
creases. The double occupancy drops even faster, re- 
sulting in an increase of 0. In the range of the transition 
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separated by a small gap A„ 



FIG. 4: Evolution of the renormalization constant Eo and 
the effective mass m/ with increasing U for e/ = — 6.0 1. The 
horizontal lines indicate the predicted transition lines to the 
(O)-TI at Eo = 1.2 1 and the (tt)-TI at E = 6.0* obtained 
from (JSJ). The solid curves are guides to the eye. 




FIG. 5: f-level occupancy, double occupancy and the measure 
of correlation O for £/ = — 6.0 i as a function of U. 



to the (O)-TI state at U ~ 1.25 1, the effect of correla- 
tions is still negligible, i.e. the transition can also be 
predicted using a mean-field decoupling of the interac- 
tion Hamiltonian Hu^. In the regime where U > 8.0 1, 
as the system approaches half-filling, correlation effects 
become very pronounced, as is indicated by the measure 
of correlation 9 taking values close to 1. Thus we have 
well-defined localized moments on the f-orbitals, and a 
small band gap A g <C t separating the weakly disper- 
sing heavy bands. This becomes manifest also in the 
spectral function A(k,u>) = —ir~ 1 ImTr G(fc, oj + iS). In 
the small frequency limit, we can expand the self-energy 



E(w) 



0) 



- uj- 



duj 



lj=0 



. This way we obtain the 



low-energy spectrum, as shown in Fig. [6^ for the weakly 
correlated (O)-TI, and in Fig. with the strongly renor- 
malized, almost flat bands of the Kondo insulating state, 




(0,ti) (0,0) 
k 



FIG. 6: Low energy part of the spectral function for £/ = 
— 6.0 1, with an artificial broadening S = 0.01. The interaction 
strength is a) U = 5.0 1, with a resulting band gap A 9 as 0.63 1 
and b) U = 8.2 1 with A 9 w 0.14 1. Clearly visible is the 
formation of the heavy bands in the strongly correlated case. 



VII. DISCUSSION & CONCLUSIONS 

We have described and studied within the DMFT ap- 
proximation a model relevant for the understanding of 
topological Kondo insulators. The concept of the topo- 
logical Hamiltonian ([3]) provides a phase diagram, which 
serves as a roadmap for the interpretation of the numer- 
ical results, as long as an adiabatic connection to a non- 
interacting state exists. The low-temperature phase di- 
agram was found to be very rich. In particular, it is 
possible to study a transition via an intermediate semi- 
metallic state between two distinct topological states, 
which is driven by the interaction, and which was not ob- 
served previously. With open boundaries, the two states, 
the (0)-TKI and the (7r)-TKI, exhibit zero-energy edge 
modes crossing at k x = and k x = ±7r respectively. 
The different nature of the two topological states was 
evidenced by establishing the adiabatic connection to the 
respective non-interacting states. Starting from a trivial 
insulating state at U = 0, we can equally observe an 
interaction-driven transition to the (0)-TKI state. This 
aspect is similar to investigations of the BHZ-Hubbard 
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model. There, the same correlation-driven transition 
from a trivial to topological insulator was foun d 14 ' 35 . 

All observed zero-temperature phases adiabatically 
connect to non-interacting states. It is however impor- 
tant to note that the topological state survives well into 
the local moment regime. This regime is characterized by 
strongly rcnormalizcd bands and a small band gap char- 
acteristic of Kondo insulators. The realization of this 
state of matter in a simple toy model allows the detailed 
study of the temperature evolution of the single particle 
spectral function on various topologies with and without 
edges. In particular such studies will provide a detailed 
understanding of the emergence of edge states as the 
temperature crosses the relevant energy scales of Kondo 
physics. Our study is not limited to the two dimensional 
case, and can likewise provide insight into three dimen- 



sional topological Kondo insulators. Again of particu- 
lar interest is the temperature dependence of the single 
particle spectral function as well as the nature, weak or 
strong, of the topological state. Studies along these lines 
are presently under progress. 
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